CIRCUIT DESIGN SUPPORT METHOD, IMPLEMENT THEREOF, AND 
CIRCUIT DESIGN SUPPORT PROGRAM 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

The present invention relates to a circuit design support method, a 
implement thereof, and a circuit design support program, and in particular, to a 
circuit design support method, a implement thereof, and a circuit design support 
program capable of automatically designing a circuit interactively through a 
10 designer and a computer without providing a circuit equation comprised of a 
differential equation. 

Description of the Related Art 
The information society of recent years has been under significant 
15 development, and increase in frequencies used for an electronic circuit has 
become increasingly prominent, for example, a computer CPU requires 
operating (clock) frequencies of the order of GHz, the use of frequencies of the 
order of GHz has also commenced in wireless communication, and the like. 
Accordingly, the development of an analog circuit has become extremely 
20 important, such as a power supply circuit for supplying power to a CPU and the 
like, and an amplifier used in wireless communication. 

The design of an analog circuit connected to an electronic circuit of which 
the operating frequencies increase, such as a power supply circuit or an 
amplifier, is categorized into the following three stages: 
25 (1 ) Decision of circuit configuration, 

(2) Decision of a circuit element value, and 

(3) Evaluation of circuit operation. 
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As a tool for supporting such design on a computer, a circuit simulator such 
as "SPICE" has been developed, which has been widely used as an important 
support tool for (3) above. However, in order to simulate circuit operation, a 
circuit element value needs to be provided, and therefore stage (2) above is 
5 essential which derives such an element value for providing a desired operation. 

A developer of an analog circuit first performs circuit analysis through 
manual calculation and then use the result of that analysis to decide the element 
value of the circuit. The circuit analysis involves many approximations, and 
therefore the derived element value includes errors. Accordingly, the developer 

10 inputs the derived element value into a simulator, and then repeatedly performs 
the above works (2) and (3) to thereby empirically adjust the element value for 
the circuit on the basis of the simulation result. Moreover, adjustment of the 
element value is largely dependent on a designer's experience, and therefore 
training the designer is also necessary. 

15 On the other hand, even if the element value is obtained with the Monte 

Carlo method using a simulator, there exists the problem that calculation time 
increases exponentially as the number of conditions to be satisfied increases. 
That is, stage (2) above that bridges the above (1) and (3) is a bottleneck in the 
development of a circuit in terms of time. For this reason, a design support tool 

20 with which even a less experienced designer can easily and rapidly design a 
circuit is strongly required. 

In consideration of this, the present inventors previously proposed a 
method for calculating an element value, in which the element value was 
numerically derived with a computer, as described in the following non-patent 

25 literature. The prior method for calculating an element value will hereinafter be 
described. 

For example, consider a circuit described with the differential equation 
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represented by the following eq.(1 ). 
(1) 



dx 



where t R, x R n , and R m represents time, an n-dimensional state space, 
5 and m-dimensional system parameters, respectively. For simplicity, we assume 
that 
(2) 

f : Rx R n x R m -> R n (t,x,\) k+ /(t,a,A) 

represents a CT-class map having periodicity of period tj in a time domain. That 
10 is, 
(3) 

/(* + t r »a:,A)~/(*,ar t A). 

represents the assumption. Furthermore, the above eq.(1) is assumed to have 
a solution of x(t) = (t, xo, ), if <t< , with respect to any initial value xo 
15 R n and all the system parameters R m :x(0)= (0, x 0 , )= x 0 . 

The assumption of the periodicity described in eq.(3) enables a C"-class 
discrete map T which returns from the state space R n to itself to be defined as: 
(4) 

T:R n ~* R n x 0 *-* T{x 0 ,\)=<i>{tT,xo,\) 
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This map is referred to as the Poincare map. 

In the case where the solution x(t) = (t, p 0 , ) of differential equation (1 ) 
has periodicity of period t T , 
(5) 



^(Po,A) =po- 
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the point p 0 R n satisfying eq. (5) is called a fixed point with respect to map T. 

In order to design an electronic circuit, constraint conditions are often 
considered. Given that the number of conditions is N (< m), by expressing each 
condition at time t = t ck as g k , the constraint condition can be expressed as: 
(6) 



0, 



fll(*cl»2CO» A) 
9N(tcN,XQ, A) 



In this case, N parameters can be prepared as the design parameter. 
Consequently, the remaining (m - N) parameters need to be provided as the 
design specification. 

As described above, the design of an electronic circuit comes down to 
solution methods of the above eqs. (5) and (6), which are algebra equations. 
These algebra equations are now rewritten as: 
(7) 



4 



F(ao,A) 



T(wo, A) - x 0 
G<ar 0 ,A) 



0, € K 



where T(x 0 , ), x 0 , and can be expressed as T(x 0 , ) = [T-i(x 0 , ), T 2 (x 0 , 

), -,T n (x 0 , )] T , x 0 = x(0) =[x 1 (0), x 2 (0), -, x n (0)] T , and =[ 1, 2, -;, 
m ] T , respectively. Further, u R N is defined as: 
5 (8) 

A„ ~ {A^i , A U 2, ■ • » , A UJ v | A u jt (k — 1 , 2, ■ ■ ■ , N) 
are unknown design parameters in A.} 

To solve design equation (7), the Newton method that is the most general 
algorithm for solving an algebra equation is used. The unknown in eq. (7) is 
expressed as u R n+N :u=[x 0 T , U T ] T > and then an iterative calculation: 
10 (9) 

is performed until the termination condition of || k+1 u - u k || < is satisfied, 

where F' R < n+N )( n+N ) represents the Jacobi matrix of F, i.e., 

(10) 
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F'(u k ) 

dxi(Q) 



1 



dx 2 (Q) 



dT 2 



dT n (u k ) &r« 

dxi(Q) 
dgi(t cU u k ) 



dxi(0) 



dx 2 (Q) 



'dXn(0) 

8gi 



dg N (t cN ,u k ) dgx 
0x7(0) dx 2 (Q) 
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In addition, k is the number of iterations and « 1 (where we assume = 1 Cf 9 ). 
In this case, u k+1 is the solution of eq.(7). 

In order to calculate eq.(9), it is required to obtain T(u k ) in eq. (10), which 
can be obtained by applying the Runge-Kutta Method to eq. (1). Further, 
respective elements of the Jacobi matrix F'(u k ) given by eq. (10) are derived by 
solving the following linear variational equations: 
(11) 



d d<p 
dt &xq 

dtdX 



df d(j> 
dx dx 0 



with 



dtp 



t~0 



ix-d\ + d\> Wlth ixh 



0 



10 Eqs. (11) can be regarded as differential equations with respect to 5 (t, u )/5xo 
and 5 (t, u k )/5 That is, if eqs.(11) are solved from 0 to t T or t ck with the 
Runge-Kutta Method, 8 (t T , u k )/3x 0 , 8 (t ck , u k )/3x 0 , 8 (t T , u k )/<9 , and 8 (t ck , 
u k )/5 can be derived. T k is a function of , and also g k can typically be 



expressed as a function of , and therefore, from these values, T k (u )/5x 0 , 
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T k (u k )/5 , gk(t C k, u k )/dx 0 , and gk(t C k,u k )/3 can be numerically derived. 

Conventionally, the unknown coefficient u has been obtained by the 
above-described calculation, and consequently u that is the design value of the 
electronic circuit has been decided. 
5 Fig. 2 is a flowchart illustrating a procedure in which a circuit design is 

performed using the above-described prior calculation procedure. The 
procedure is comprised of 5 steps: (1) the circuit configuration decision step of 
deciding electronic circuit configuration; (2) the circuit equation formulation step 
in which a circuit designer manually decides circuit equation (1) on the basis of 

10 the decided circuit configuration; (3) the variational equation derivation step in 
which a circuit designer derives the variational equations (11) from the circuit 
equation; (4) the constraint condition derivation step of deriving the constraint 
condition (eq. 7); and (5) the step of deriving a circuit element value and a 
waveform by solving conditional equation (7) with circuit equation (1) and 

15 variational equations (11 ), and if a satisfactory solution is not obtained in step 5, 
then an operation returns to step 1 . 

Non-patent literature 1: IEEE TRANSACTION ON CIRCUITS AND 
SYSTEM-1: fundamental theory and applications, vol. 49, No. 7, July 2002. 

As described above, the prior method is configured to derive a voltage or 

20 current waveform x(t) = (t, xo, ) by providing the circuit equation in an explicit 
form following the consideration of circuit configuration and then solving the 
circuit equation. However, it is not easy to automatically derive the circuit 
equation with a computer, and therefore, there exists a problem of having to 
depend on manual calculation. Furthermore, every time a circuit configuration 

25 or the characteristic of a circuit element is modified, the circuit equation has to be 
reformulated, resulting in a large expenditure of time. In other words, in the 
prior method, a series of works, i.e., deriving the circuit equation from the circuit 
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configuration and then deriving the variational equations constitute barriers to an 
automated design. 

It is therefore a purpose of the present invention to provide a circuit design 
support method, a implement thereof, and a circuit design support program with 
5 which a circuit element value and a waveform can be derived automatically with 
a computer. 

SUMMARY OF THE INVENTION 
A circuit design support method, a implement thereof, and a circuit design 
10 support program according to the present invention is characterized by a circuit 
design support for calculating a circuit element value of an analog electronic 
circuit to be designed, by instructing a computer to execute a program describing 
an iterative calculation equation of the Newton method including the Jacobi 
matrix, wherein the program describing the calculation equation in which an 
15 approximate equation is substituted for a partial differentiation that is an element 
of the Jacobi matrix is used, and the element as the partial differentiation of the 
Jacobi matrix is obtained from a waveform observed in response to a circuit 
configuration of the analog electronic circuit. 

The present invention is further characterized in that a steady state of the 
20 analog electronic circuit is obtained simultaneously with the calculation of the 
circuit element value. 

The present invention is further characterized in that statistical data are 
provided as a constraint condition of the analog electronic circuit. 

This configuration does not need describing a circuit equation of an analog 
25 electronic circuit as a differential equation. It also does not need applying the 
Runge-Kutta method to the linear variational equation to calculate the element 
as an object of the above-mentioned partial differentiation. 
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The present invention eliminates the needs for obtaining the element of the 
partial differentiation of the Jacobi matrix from linear variational equations and 
also for describing the circuit equation in an explicit form, so that manual 
calculation by a designer becomes unnecessary. Furthermore, statistical data 
5 can be provided as a constraint condition of an analog electronic circuit, so that 
the present invention can meet wider design requirements and facilitate setting 
constraint conditions. 

DETAILED DESCRIPTION OF THE INVENTION 
1 0 One embodiment of the present invention will hereinafter be described with 

reference to Fig. 1 and the equations used for describing the prior method. 

Fig. 1 is a flowchart illustrating a procedure of a circuit design support 
method according to one embodiment of the present invention. The circuit 
design support method is comprised of step 11 of deciding a circuit configuration, 
15 step 12 of setting conditional equation (7), and step 13 of inputting the circuit 
configuration decided in step 11 into a circuit design support implement 
(hereinafter referred to as a "circuit simulator") to derive a waveform as 
described later and of solving the conditional equation input into the circuit 
simulator to thereby derive a circuit element value (a design value). In step 13, 
20 eq. (1 ) is also solved at the same time, and therefore a waveform equation when 
a condition is satisfied, i.e., a steady state of an analog electronic circuit is also 
derived. 

Next, a circuit element value design support method according to this 
embodiment is described in comparison to the prior method. 
25 In this embodiment, differential equation (1) is not set in an explicit form 

unlike the conventional method. That is, in a certain deterministic system (a 
system that can be described with differential equations) with respect to a 
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variable x, it is assumed that a solution (output waveform) of x(t) = (t, x 0 , ) 
can be observed with respect to any initial value x 0 R n and all system 
parameters R m :x(0)= (0, x 0 , )= x 0 , where x(t) has periodicity of period t T 
in a time domain. 

5 Even under this assumption, following a consideration similar to the prior 

one causes the design of an electronic circuit to come down to the solution 
methods of eqs. (5) and (6), and rewriting these eqs. allows eq. (7) to be 
obtained and eq. (8) to be defined, where, in this embodiment, eq. (6) is not 
necessarily a function of time t C i ~ t Cn , but functions of x 0 and . Consequently, 

10 in the prior method, only a condition for a time response at a certain time t c 
derived from the circuit equation can be included in a constraint condition of an 
analog electronic circuit, whereas, in this embodiment, a waveform observed on 
a circuit simulator, and statistical data such as the maximum value, the minimum 
value, and the mean value of observations calculated from the waveform can be 

15 used as constraint conditions. For example, statistical values of time 
responses . derived from the circuit equation, observations in a domain other 
than a time domain, such as the Fourier domain of , and statistical values of 
the observations can be used as the constraint conditions. 

As a solution method of eq. (7), this embodiment also employs the Newton 

20 method, and eq. (9), which is an iterative calculation, is iteratively calculated until 
the termination condition is satisfied. In this calculation, the use of the Jacobi 
matrix (10) is also the same as in the prior method, and u k+1 is the solution of 
eq.(7). 

Then, in this embodiment, T(u k ) for calculating eq. (9) is derived from 
25 observed x(t) = (t, x 0 , ). On the other hand, the elements of the Jacobi 
matrix F'(u k ) (10) cannot be solved with variational equations (11) since circuit 
equation (1) is not described in an explicit form in this embodiment, differently 
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from the prior method. 

Therefore, in this embodiment, instead of the variational equations, a 
partial differentiation that is each of the elements of the Jacobi matrix is replaced 
by an approximate equation based on a definition. That is, when 
5 (12) 

is defined, the following approximate equation (13) is substituted for a value of 
the partial differentiation that is each of the elements of the Jacobi matrix. 
(13) 

_ - g(u k ) 

dui € 

10 

where is an infinitesimal coefficient g(u j) can be derived by substituting u i 

for a parameter u of g(u) in the same deterministic system (circuit configuration) 

as that for deriving g(u). 

From the above calculations, the unknown coefficient u can be obtained, 
15 and consequently the design value u is decided, whereby the design of the 

electronic circuit becomes possible. 

According to this embodiment described above, since a partial 

differentiation is derived only from an output waveform, a variational equation is 

not necessary, and also, since a circuit equation is not given in an explicit form 
20 but in an implicit form, a manual circuit design process becomes unnecessary 

and performing all the circuit design processes becomes possible interactively 

with a computer. 

Furthermore, even in the case of design change of a circuit element, i.e., 
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even in the case where a deterministic system is changed, a parameter for 
obtaining a periodic orbit can be easily derived in each case and the change can 
be easily responded to since an approximate value is substituted for a value of a 
partial differentiation. Still further, in comparison to Fig. 2 illustrating the prior 
5 method, this embodiment is configured to eliminate the formulation of the circuit 
equation (step 2) and the derivation of variational equations (step 3), and to 
substitute an approximate equation for a value of a partial differentiation, so that 
a large change in an internal configuration of a simulator is not required. 
Consequently, a circuit simulator in this embodiment can effectively utilize 
1 0 properties of other existing simulators such as "SPICE". 



Industrial applicability 
The present invention enables a computer to calculate a circuit element 
value without providing a circuit equation when an electronic circuit is designed, 
15 so that the circuit design support implement enabling an electronic circuit to be 
designed in a fully automatic mode can be configured, which is useful for 
designing an analog electronic circuit and an mixed analog/digital circuit, 
especially a switching circuit. 



20 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a flowchart illustrating a procedure of a circuit design support 

method according to one embodiment of the present invention. 

Fig. 2 is a flowchart illustrating a procedure of a prior circuit design support 

method. 

25 

Explanations of letters or numerals 
11 : Decision step of circuit configuration, 
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12: Setting step of constraint condition, and 

13: Derivation step of design value and waveform. 
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